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An explanation of the solar transition region 

Philip Judge 

High Altitude Observatory, National Center for Atmospheric Research^, P.O. Box 3000, 

Boulder CO 80307-3000, USA 

ABSTRACT 

Prompted by high resolution observations, I propose an explanation for the 
40+ year old problem of structure and energy balance in the solar transition 
region. The ingredients are simply cross-field diffusion of neutral atoms from 
cool threads extending into the corona, and the subsequent excitation, radiation 
and ionization of these atoms via electron impact. The processes occur when- 
ever chromospheric plasma is adjacent to coronal plasma, and are efficient even 
when ion gyro-frequencies exceed collision frequencies. Cool threads - fibrils and 
spicules perhaps - grow slowly in thickness as a neutral, ionizing front expands 
across the magnetic field into coronal plasma. Radiative intensities estimated 
for H La are within an order of magnitude of those observed, with no ad-hoc 
parameters - only thermal parameters and geometric considerations are needed. 
I speculate that the subsequent dynamics of the diffused material might also 
explain observed properties of trace elements. 

Subject headings: Sun: atmosphere - Sun: chromosphere - Sun: transition region 
- Sun: corona - Sun: magnetic fields 
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Introduction 



The upper transition region (henceforth, "TR" ) - plasma with electron temperatures in 
the range 2 x 10^ ^ ^ 10^ K, is adequately described by field-aligned thermal conduction 
down from t he corona. The lower TR (10^ K < Te < 2 x 10^ K) however, is not so easily 



understood flCabrielll 19761 : 1 Jordanlll980l ) . Models dominated by field-aligned heat conduction 



produce too little emission from the lower TR by orders of magnitude, a problem already 



rected conductive flux of F, 



cond 



XQ6 gj^gg 2 g 



[e.g. 



Jordan 


1980; 


Athav 


1981) 



The National Center for Atmospheric Research is sponsored by the National Science Foundation 



-2- 



et al. (2002 and earlier papers in the series), henceforth "FAL" showed that energy balance 
can be achieved through field-aligned (ID) diffusion of neutral hydrogen and helium atoms. 
The neutral atoms diffuse into hot regions, radiate away much of the coronal energy, and 
can reproduce the H and He line intensities. 

The problem might be considered by some as solved, in principle. But there exists 
the serious and nagging proble m of the pecul iar spatial relationship between the observed 
corona, TR and chromosphere (jFeldmanlll983l ) . Feldman and colleagues have since analyzed 
man y observations, conclu ding that the lower TR is the rmally disconnected fr om the corona 
(e.g. iFeldman et a/.ll200ll . and references therein). Yet iFontenla et al\ ( Il990l ) declared that 



"The above [i.e. their] scenario explains why (as noted by Feldman 1983) the structure of 
the transition region is not clearly related to the structures in the corona" . That the debate 
still rages is evidenced by advocates for "cool loop" models in which lower TR radiation 
originates from loops nev er reaching coronal temperatures and having negligible conduction 
( jPatsourakos et a/.l 120071 . and references therein, henceforth "PGV"). Here I propose a dif- 
ferent scenario, prompted by new data and analyses which show that neither cool loops nor 
field-aligned processes adequately describe the La chromospheric network. I speculate that 
other TR lines might also be accounted for. 



2. A new scenario 



La network emission, at 0'.'3 resolution appears mostly as threads of relatively uniform 
intensity, of 5-10 Mm length and IMm diameter (PGV). PGV argued that "the different 
appearance the TR has in the quiet Sun [i.e. network] is suggesting that the bulk of its 
emission comes from structures other than the footpoints of hot loops". Convolved La 
images from P GV appear to cor r espond to those see n in rn any other TR lines at lower 



resolution (e.g. ICurdt et a/.ll200ll ). IJudge and Centend (120081 ) showed, using magnetic field 



measurements from Kitt Peak, that much of the network La emission originates in long 
spicule-like structures lying along the lowest few Mm of magnetic field lines extending into 
the corona, but that plage emission may correspond to the thin footpoints as suspected by 
PGV and modeled by FAL. Even in plages, on sub-arcsecond scales, field- aligned threads 
of cool plasma (fibrils, spicules), extend into the low corori a forming non-planar thermal 
interfaces between hot and cool plasma (IBerger et a/.lll999l ). Prompted by these data, I 
examine the dif fusion of neutral particles into the corona, across magnetic fields (following 
a suggestion by lPietarila and Judge! l2004l ). 



Consider a straight cylinder of cool, partially ionized material embedded in a hot corona, 
of radius r^. Length l^. ^ of the tube contains cool plasma in contact with the hot corona. 
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The magnetic tube is of length L ^ Ic, mostly containing coronal plasma. Tube parameters 
are given in Table [TJ Note that the neutral density greatly exceeds other densities. The 
chosen geometry is typical of values found by PGV, and thermal parameters are typical of 
the quiet Sur0. The plasma is assumed to be in a low plasma-/? regime. 



2.1. Initial diffusion, relaxation, radiation 



Imagine an injection of dense neutral material into the tube footpoint by some chro- 
mospheric process. The tube surface acts as a semi-permeable membrane. Neutral particles 
travel freely between collisions, but ions gyrate about magnetic field lines with gyro radii 
orders of magnitude smaller than mean free paths ("mfp"s, Table [2]). Ions and electrons are 
essentially frozen to field lines, but neutrals can diffuse across field lines almost as efficiently 
as along them, and find themselves impacted by hot electrons and protons. 

Table [2] lists time scales for kinetic pro cesses for a "c ool" hydrogen atom e mbedded in 
a hot corona of T = lO^K, usin g data from lAUenl (jl973l ). iHansteen et al\ ( 119971 . henceforth 



HLH), and [Gilbert et al\ ( 120021 ). A hydrogen atom crossing the boundary encounters other 
diffusing hydrogen atoms and hot protons and electrons. Statistically, the first interaction is 
a collision with a coronal proton, involving the exchange of energy and (~ 50% of the time) 
an electron (charge transfer, "CT"). Charge transfer y ields an exchange of momentum (180° 
change in direction) but little exchange of energy (e.g. IOsterbrocklll96ll ). The kinetic energy 
exchanged is ~ ^kT^, shared between them after two such collisions (I use subscript "/i" to 
denote hot and "c" cool plasma). The CT cross section is roughly independent of energy, 
so the "warm" neutral atom has a ~ 1 — = 0.63 probability of staying within the hot 
plasma. Assuming that it does so, after ti^. ^ 8tct s it becomes ionized by impact with a 
hot electron. Once free, the electron will not readily recombine with a proton (time scale 
Tki ~ 3 X 10^ s). At Te = Th = lO^K, the time tu needed for electron impact excitation 
of the n = 2 levels and (rapid) emission of a La photon is comparable to ri^. Thus, of the 
ionized neutral atoms, ~ 50% will have emitted a La photon, the energy supplied by coronal 
electrons and protons. Because there are relatively few hot particles, their thermal energy 
limits the number of neutral atom ionizations and excitations. 



From kinetic theory, the flux density of neutral hydrogen atoms initially crossing the 
boundary into the corona is ^nJJc ~ 2 x lO"*^^ particles s~^ cm~^, where v = \/%kTcli\m. 



The kinetic energy per "hot" proton is \kTh which is shared roughly equally after two CT 



do not adopt the higher temperatures of cool loops used by PGV, because here the corona and dynamics 
supplies all the energy for La emission. 
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collisions with a neutral, producing a population of rih <^ Uc "warm" neutrals with T ~ Th/2. 
After a few, say m more collisions (time rriTcT later), all the proton energies are the larger 
of ~ and |fcTc, and mnh of the Uc neutrals have suffered a proton impact. (A time 

of ric/nh = 100 times tct ~ Is is required before all neutrals have been impacted). The 
warm neutrals relax via collisions with the cool neutrals. 

The initial electron evolution is largely determined by inelastic collisions with hydrogen: 
each hot electron typically has sufficient energy to excite and ionize 5 neutral hydrogen atoms, 
which takes ~ Tri^ ~ 0.6s. (Electron-electron collision times are < 0.05 s). The electrons 
lose energy e = 5nh{I + E)e per unit volume at the rate 

- ^ — ~ 0.13 erg cm"^ s"\ (1) 

t YTiK 

leading to a cooling time of < 0.4 s. (A lower limit applies because tails of the Maxwellian 
distributions can increase excitation/ionization rates). Of this energy a fraction = | is 
emitted in La. The radial flux density of La radiation from this neutral "sheath" is 

/ > ^^A ^ ^ev^^^f erg cm'^ s-\ (2) 

where A is the sheath thickness at time t (60rc"r); and ti^*-^^ = A/t is the diffusion speed. 
For a random walk, Ac ~ ^^/QOtctVc ~ 3.3 x 10^ cm, for warmed neutrals A^ ^ 1.8 x 10^ 
cm {v^^^'f = 0.57 and w^*-^-^ = 3.2 km/s respectively; the factor | accounts for the random 
direction of the "walk"). As a rough estimate, I take u*-^-^ Sv^^'^^ km/s: 

/ ^ ^eSvf^^ ^ 5.6 X 10^ erg cm'^ g-i (^3) 

Thus, / is initially just a fixed fraction of the local coronal energy density multiplied by the 
diffusion speed. The specific intensity / equals //vr when the line is optically thick and all 
the radiation scatters away from the solar surface. (Photon mfps for La in the sheath are just 
10^ cm). This estimate of / is a factor of 100-300 below measured values of (1.8 — 5.6) x 10^ 
ergs cm~^ s~^ sr~^ in active n etwork threads (PGV) and 30 below average network intensities 



( iVernazza and Reeved 1 19781 ). But, as will be made clear below, this is an under-estimate. 
Similar estimates of intensities for H L/3 and the 584A line of He I, relative to La, are quite 
reasonable, recognizing that L(3 is optically thick across the sheath. 



2.2. A multi fluid calculation 



To examine the evolution at later times, multi-fluid equations for conservation of mass, 
momentum and energy were solved as functions of time and distance x across the fleld lines 
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following ISchunkl (119771 ) and HLH. Just electrons, protons and neutral hydrogen atoms were 
treated. Cartesian geometry is used because the diffusion region is much thinner than the 
tube. I assume that electrons are strongly tied to protons, so that their densities and ffuid 



velocities are equal 



Tip, tie 



Ur 



charge and electrical currents are neglected). The 



conservation equations used for mass, momentum and energy density for the fluid of species 
s are 



dn^u. 



dt 



dus d 6ns 
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+ —{msUsul+ps + rff } + F 



dEs 
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+ Q 
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(4) 
(5) 
(6) 



dt dx^ ' ' "^'^ ' ^ 5t 
No conservation equation is used for the heat flux since here it is treated as df using the 
mfp approximation. Above, F is a body force term (gravity, Lorentz force for example). 
Eg = \nskTs + ^nisnsul, Ps = nghTg, and the |^ are non-linear collisional terms. Q and L 
are the energy gains and losses respectively, where I adopt Q = 1.67 x 10~^^?T,e^ipe~"^^/^°°° 
erg cm^ s~^ to maintain a chromosphere against losses L (HLH), and L includes latent heat 
and La radiative losses computed explicitly from the collisional terms. 

The diffusion terms dg (not included by HLH, except for the heat flux) require care espe- 
cially for the dynamics of the proton fluid. For individual protons and electrons, the momen- 
tum equations are dominated by the Lorentz force. Their cross-field motion on timescales 
short compared with collision times is circular with frequency Ug = CgB /mg. On longer time 
scales the summed (fluid parcel) momenta can change only after a collision. The net effect of 
the Lorentz force is thus to limit the cross-field displacement of charged particles to a single 
gyro radius Tg = Vg/cOg in collision time Tg instead of the collisional mean free path = VgTg. 
Thus, a simple recipe for calculating cross-field transport via the fluid equations is to set 
both F and ^ terms to zero in the proton momentum equation, and modify the dp terms to 
account for the red uced displacemen ts. Field-free diffusion is described by equations (4.41), 
(4.46) and (4.52) of lOombosil (|l994j ): 



1 d _ 

3 dx^ ^ 
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M 



-As— {mgngUgVs} , 



dt 
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12 'dx ^ 



7tt3 



}, (7) 



For charged particles, A, mus t be replaced by A* = As/(1 + ujgTg)'^ (following the above 
argument, see lBraginskiilll965l . eqs. 4.37, 4.40). Note that, written in terms of Tg, df yields 
the widely used "Spitzer" thermal conductivity parallel to the field, and the ion-dominated 
conductivity perpendicular to the field. The net effect for ujgTg ^ 1 is that only the neutral 
fluid diffuses efficiently across the field- the charged fluid evolves mostly via the collisional 
coupling to the neutrals (via 6Mg/6t), and to a lesser degree to the small dg terms. 
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The variables (n^, Us, Eg) for electrons, protons and neutral H atoms, functions of [t; x), 
were initialized according to table [TJ Only 7 variables were solved since it is assumed that 
= np and = Up. The equ ations were integrated in tim e using MacCormack's method 



to include the coUisional terms (iGriffiths and Highanj|l999l ). For the first three points near 



X = 0, the variables were held fixed to their initial values, maintaining the same chromo- 
spheric conditions there. Figured] shows conditions several seconds after the beginning of the 
diffusion process. Pressure gradients drive neutrals into the corona against friction forces, 
thus the diffusion speed, measured by tracking the steep temperature rise, is ~ 0.8 kms~^, 
far below the thermal speed. The computed flux density of La is ~ 5 x 10'^ ergscm"^ s~^, 
and is roughly constant in time. It is some 10 x higher than the simple kinetic result above, 
because of the nonlinear dynamics: (1) the densities become higher in the corona, (2) flow 
energy is converted to heat, (3) the La losses/latent heat ratio is higher (the photons are 
created at electron temperatures lower than the initial coronal temperature). / is computed 
to be just a factor of 10-30 below observed active network thread intensities, and 3 below 
average network intensities. 

A calculation with twice the coronal density, more appropriate for active network, yields 
smaller diffusion speeds and La fluxes which are just 1.7 times higher. EUV/X-ray coronal 
intensities scale with (density)^, and so would be a factor of four brighter. This non-linear 
relationship is an important property of the calculations. 



3. Discussion, speculations 

Based upon observations of spicules and other fine, thread-like structures on the solar 
disk, it is clear that non-planar thermal interfaces exist at the base of the corona, and that 
the morphology of the TR emission from such interfaces cannot be explained by field-aligned 



particle transport at the base of coronal loops, in contrast to the claims by iFontenla et al 



(Il990l ). The picture proposed here uses unspecified mechanisms in the chromosphere to 
maintain a reservoir of cool mostly neutral plasma directly adjacent to hot coronal plasma. 
The cylindrical geometry, inspired by observations, presents a large surface area (per unit 
volume) of contact between cool and hot plasma. The chromosphere supplies mass via neutral 
diffusion across the surface to a thermal boundary layer, and the corona supplies energy to 
the neutral particles. The originally neutral particles drain energy from the corona by latent 
heat of ionization and by inelastic collisions leading to strong La emission. The diffusing 
layer propagates outwards, emitting radiation like the boundary of a wild fir^, into the 



^Secchi in 1877 described the chromosphere as a "burning prairie", but in a different sense. 
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corona until either the supply of neutral mass or coronal energy dries up. The present 



proposal is related to models invoking cross- field heat conduction (IRabin and Moorelll984j : 
Athaylll990l ). This effect is included here (via df), but it is far less efficient at moving heat 



to cool plasma than diffusion is at moving neutral atoms to the coronal heat. 

The calculations presented here fall short of accounting for the large radiative flux of La, 
by factors of ~ 10. However, the calculations miss important additional sources of energy in 
the corona: thermal and gravitational potential energy. The cool threads extend only a few 
Mm into the corona, and form just the lower parts of a much larger coronal structure. The 
diffused cool material is thus subject to parallel transport (heat conduction, diffusion) which 
will transfer heat from the overlying coronal plasma to the diffused material. Spicules formed 
by ejection from the chromosphere will have their entire length exposed to this energy flux, 
because the lowest parts of the spicules diffuse flrst into the corona- the diffusion fronts are 
not exactly parallel to fleld lines. Coronal plasma along connected field lines contains L^nkT 
erg cm~^, where L is the pressure scale height (~ 50Mm) or loop length. Since L ^ l^. the 
energy available for La radiation would be L//c > lOx larger than computed above, more if 
the tubes expand with height. I speculate that cross-field diffusion and subsequent parallel 
conduction might bring theoretical and observed intensities values into agreement. The time 
needed to conduct this energy must lie between the electron sound speed Ce as L/ce ~ 13 s, 
and ~ 10'^ s, an upper limit obtained from the thermal energy divided by the conductive flux 
for a uniform temperature gradient. Gravitational potential energy might contribute to the 
heating and dynamics of the sheath as the diffused material cools the corona and adds mass, 
such that vertical pressure balance no longer is expected. It may be that larger red-shifts 
would be expected where magnetic flelds are more vertical, i.e. directly ov er the magnetic 



network. This expectation is not in disagreement with results found by iMcIntosh et al. 



(120071 ). However, little more can be said without solving the 2D multi-fluid conservation 
equations including parallel heat conduction and cross-fleld diffusion, beyond the scope of 
this letter. Such calculations will also show if the emission lines of trace species (ions of 
carbon, oxygen etc. in the TR) can be explained. 

Cool threads are observed in different coronal environments (PGV)- their intensities 
appear to vary relatively little compared with the embedding coronal intensities. This fact 
is part of Feldman's (1983) claim that TR emission is energetically disconnected from the 
corona. The calculations presented here indeed produce a non-linear relationship between 
La and coronal brightness. The La intensities scale with the local coronal energy density 
and with the diffusion speed. But the EUV and X-ray radiation emitted by the corona itself 
vary with density^ and peaked functions of temperature along lines of sight different from 
the direction of fleld lines into the sheath. The scenario might therefore explain most of 
the observed puzzling facets noted by Feldman and colleagues, yet still maintain a strong 
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energetic link between the c orona and TR, and t her eby resolve a long-st anding debate (see 
the different perspectives of iFeldman et a/.ll200ll and lWikst0l et a/.lll998l . for example). 

To see if the scenario survives scrutiny, more observations of chromospheric fine structure 
and its relation with the corona and TR would be as important as numerical modeling work. 

I am grateful to Tom Holzer, Scott Mcintosh and the referee for comments. 
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Table 1. Tube properties 



Quantity 


Inside 


Outside 




(cool) 


(hot) 


radius cm 


5 X 10^ 




length Ic cm 


5 X 10^ 


> 10/c 


T K 


10^ 


10^ 


uh cm~^ 


8 X 10^° 





Up cm~^ 




4 X 10^ 






4 X 10^ 


p dyncm"^ 


0.11 


0.11 


Magnetic field strength B G 


10 


10 


B'^/Stt dyncm"^ 


3.8 


3.8 
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Table 2. Plasma conditions 



Quantity Units 



scaling 



notes 



Initial corona 
Th 

rih 

P 

B 

Up 
fgyro 



^p'^pp 



K 

cm" 
cm" 
cm" 
G 

s-i 
km 
s 



10^ 

8.0 X 10^ 

4.0 X 10^ 

1.1 X lO'i 
10 

2.8 X 10-2 
9.6 X 10^ 
1.5 X 10-^ 
1.6 

1.5 X 10^ 
5.0 X 10-2 

8.0 X 10^ 
13 

Tnn S 1.4 X 10^^ 

hot protons impacting hydrogen atoms 

Tpn{CT) S 1.0 X 10-2 

H atom mfp km 6.5 x 10"^ 
cool hydrogen atoms impacting protons 



chromospheric tube 
Tc K 
V km s" 



rir 



cm 



np 



(CT) 



proton mfp 



s 

km 



^pTnp 



8.0 X lO--' 
5.8 X 10-3 
7.7 

hot electrons impacting H atoms 

9.5 X 10-2 
8.2 X 10-2 
4.0 X 10^ 



Tl2 
Tlk 
Tkl 



rjpl/2 



n~ 



12.-1/2 
12.-1/2 



ly-1/2 



"CT" — charge transfer 



n^-lTe~'/'el3.6e/fcTe 



12^+1/2 



excitation of n = 2 level 
ionization 

radiative recombination 
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Note. — Tab refers to the time taken for a particle of type b to be impacted by a sea of 
particles of type a, except where noted. 
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Fig. 1. — Conditions a few seconds after hydrogen is allowed to diffuse across field lines into 
coronal plasma. The abscissa is distance x across the field lines, the initial cool flux tube 
extends from zero to the dashed line. 



